Mechanical and chemical spinodal instabilities in finite quantum systems 
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Self consistent quantum approaches are used to study the 
instabilities of finite nuclear systems. The frequencies of mul- 
tipole density fluctuations are determined as a function of 
dilution and temperature, for several isotopes. The spinodal 
region of the phase diagrams is determined and it appears 
that instabilities are reduced by finite size effects. The role of 
surface and volume instabilities is discussed. It is indicated 
that the important chemical effects associated with mechani- 
cal disruption may lead to isospin fractionation. 
PACS numbers: 21.65.+f, 25.70.Pq, 21.60.Ev 

The dynamics of first order phase transitions is often 
induced by instabilities against fluctuations of the order 
parameter. For instance, mechanical instabilities lead to 
liquid-gas phase transitions and chemical instabilities in- 
duce spinodal decomposition of binary alloys. In violent 
heavy ion collisions, nuclear matter may be quenched in 
the coexistence region of the nuclear liquid-gas phase dia- 
gram. Then, the observed abundant fragment formation 
may take place through a rapid amplification of spinodal 
instabilities. New experimental results pleading in fa- 
vor of such a spinodal decomposition have recently been 
reported From the theoretical point of view, the 

spinodal instabilities in finite systems have been mainly 
studied within semi-classical or hydrodynamical frame- 
work However, since the relevant temperatures are 
comparable to the shell spacing and the wave numbers 
of the unstable modes are of the order of Fermi momen- 
tum, quantum effects are expected to be important, as 
stressed in ||[)). 

In this work, we present a fully quant al investigation 
of the spinodal instabilities and related phase diagrams 
of finite nuclear systems. We determine frequencies and 
form factors of the unstable collective modes of an ex- 
cited expanded system by linearizing the time dependent 
Hartree Fock (TDHF) equation. We carry out applica- 
tions for Ca and Sn isotopes. The quantum nature of the 
system is responsible for many dominant features such 
as the fact that the first mode to become unstable is the 
low lying octupole vibration. Slightly diluted systems are 
unstable against low multipole deformations of the sur- 
face, which may be associated with asymmetric binary 
or ternary fission processes. We extend our analysis to 
charge asymmetric systems and show the importance of 
chemical effects in the spinodal fragmentation process. 

In TDHF theory, the evolution of the one-body den- 



sity matrix p(t) is determined by, ihdtpit) = [h\p], p(t)], 
where h[p] = p 2 /2m + Z7[p] is the mean- field Hamiltonian 
with U[p] as the self-consistent potential. To investigate 
instabilities encountered during the evolution of an ex- 
panding system, one should study the dynamics of the 
small deviations Sp(t) around the TDHF trajectory p(t) 
]To[ . It is more convenient to carry out such an investi- 
gation in the "co- moving frame" of the system, which is 
described by the density matrix, p(t) = (t) p(t) U (t) 
where U (t) = exp[-jr\tQ] with Q as a one-body opera- 
tor. For example, in our case Q could be a suitable con- 
straining operator for preparing the system at low den- 
sities and A is the associated Lagrange multiplier. Then, 
the TDHF equation in the moving frame transforms into, 



ifi^p(t) = [h(t)-\Q, P (t)\, 



(1) 



where the mean-field Hamiltonian in the moving frame 
is given by, h(t) = U 1 " (t) h[p] U (t) . The small density 
fluctuations Sp(t) in the moving frame are determined by 
the time-dependent RPA equations, 



ihpp = \h(t) 



XQ,Sp] + [SU,p}. 



(2) 



Here, we consider the early evolution of instabilities 
in the vicinity of an initial state po determined by a 
constrained Hartree- Fock solution [h(0) — AQ, po] — 0, 
where h(0) = h[po] = ho is the mean- field Hamiltonian 
at the initial state. Then, small density fluctuations are 
characterized by the RPA modes p v and the associated 
frequencies uj v . Incorporating the representation \i >, 
which diagonalizes ho — XQ and po with eigenvalues t{ 
and occupation numbers n$, the equations for the RPA 
functions < i\p u \j >= p z J become, 



hu vP l J = (a-cj) pl J +2^(n r m)Vii ikj p k v 

kl 



kl 



(3) 



where Vg^kj =< i\dU/dpik\j > denotes the residual in- 
teraction [l^Jl^]. When the frequency of a mode drops 
to zero and then becomes imaginary, the system enters 
an instability region. 

To performe an extensive study of instabilities we may 
parametrize the possible densities po either by a static 
Hartree-Fock (HF) calculations constrained by a set of 
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collective operators [12 , or using a direct parameteriza- 
tion of the density matrix. We follow the second ap- 
proach by introducing a self-similar scaling of the HF 
density as suggested by dynamical simulations. 

We solve the HF equation for the ground state 
[Hhf, Phf] = 0, leading to the single-particle wave func- 
tions \cpi > and the associated energies S{. We intro- 
duce the density matrix at a finite temperature T as 
Phf [T] = 1/ (1 + exp ((h HF -e F [T]) /T)), where e F [T] 
is the corresponding Fermi level that is tuned in or- 
der to get the correct particle number. We perform a 
scaling transformation, R [a] , which inflates the wave 
functions in radial direction by a factor a according 
to < r|i?[a]|<£ >= a -1 / 3 < r/a\ip > . We define 
the density matrix for a hot and diluted system by 
po [a,T] = R[a] phf [ a<2 T] R^ [a]. The associated con- 
strained Hamiltonian is thus ho [a] = orR [a] hnF a 
so that the constraint can be identified as —XQ a = 
ho [a] — h[po [a, T]]. By construction [ho [a] , po [a, T]] = 
0, so that they can be diagonalized simultaneously. The 
eigenstates of the constrained Hamiltonian are given 
by \i >= R[a]\(fi >, and the corresponding energies 
and occupation numbers are = €i/a 2 and rii = 
1/ (l + exp ((si-ep [ct 2 T~\^ /a 2 T)), respectively. 
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FIG. 1. Contour plots of the isoscalar strength functions 
associated with the multipolarity L = 2 — 5 as a function of 
the dilution parameter a and the collective energy of the mode 
E u — chuju (c = 1 for stable modes, —i for unstable modes) 
for 40 C a (top) and 120 Sn (bottom). The strength values are 
given in percentage of the EWSR: 2 (yellow), 15 (orange) and 
30 (red). 

We perform the HF calculations in the coordinate rep- 



resentation using the Skyrme force SLy4 [^3[ . We note 
that the isospin symmetry is already broken at the HF 
level. The particle and hole states are obtained by diag- 
onalizing the HF Hamiltonian in a large harmonic oscil- 
lator representation Q, which includes 12 major shells 
for Ca isotopes and 15 for Sn. We apply the scaling 
and heating procedures described above to the density 
matrix, and calculate the residual interaction in a self 
consistent manner. We solve the RPA eq.(^) by a direct 
diagonalization using a discrete two quasi-particle exci- 
tation representation 

The top part of figure 1 shows calculations performed 
for 40 Ca. Top panels shows contour plots of the isoscalar 
strength function associated with the isoscalar operator 



LM 



^2 i=i ' 



j Ylm, with multipolarity L = 2 — 5, as 



a function of the dilution parameter a. In the stable do- 
main, the energy associated with the dominant isoscalar 
strength decreases as dilution becomes larger, and at 
a critical dilution it drops to zero. At larger dilution, 
the system becomes unstable, and for each multipolarity, 
one or two unstable modes appear. This is illustrated 
in the bottom panel, where the "energy" of the mode 
E v = —ihuj v is plotted as a function of the dilution. 

It is seen that at low dilutions, around a = 1.2, only 
the octupole mode becomes unstable. In general, density 
fluctuations with odd multipolarity become unstable at 
relatively smaller values of dilution than those for even 
multipole fluctuations. This is a genuine quantum effect, 
due to the fact that the majority of the particles have 
to jump only one major shell to produce an odd natural- 
parity particle-hole excitation while twice this energy is 
required for an even one. In nuclei at normal density this 
makes the first 3~ a strongly collective low lying state; 
in diluted system this state is the first to turn unstable. 
The second feature is that large multipole deformations 
are hardly becoming unstable. This is due to surface 
and quantum effects that prevent the break-up of such a 
small system into several fragments. As a consequence, 
the fastest growth time, r v = hj\E v \ « 28 fm/c, occurs 
for L = 2 at the dilution a=1.8. However, deep inside 
the instability region, the octupole mode is almost as 
unstable as the quadrupole mode. 

Results of similar calculations performed for 120 Sn are 
shown in the bottom part of figure 1. Also in this case, 
the octupole mode becomes unstable first, at a ~ 1.1. 
Large multipoles are more unstable in 120 Sn than in 40 Ca, 
since the system is larger and can afford larger multipole 
deformations. This is the reason why the smallest growth 
time in 120 Sn occurs for L = 3. 

It is useful to study the behavior of the RPA solu- 
tion in the coordinate space, p+ (r) = p+ (r, r), for 
protons and for neutrons separately. Since the unstable 
modes have large isoscalar strength, protons and neu- 
trons mostly move in phase. However, because of the 
isospin symmetry breaking induced by the Coulomb force 
and by the initial asymmetry of the system considered, 
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neutrons and protons oscillations have different ampli- 
tude and shape. Figure 2(a) shows the radial dependence 
of the form factor associated with the unstable octupole 
mode, at the dilution a = 1.5, for protons (dotted line), 
neutrons (full line) and the sum (dashed line), in Sn iso- 
topes. Contour plots of neutron (proton) perturbed den- 
sities, Po p ,n( r ) + Pp,n( r )i are a ^ so snown in parts b (c) of 
figure 2. 
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FIG. 2. (a) Radial dependence of the form factor associated 
with the unstable mode with L=3, at the dilution a — 1.5, for 
protons (dotted line) , neutrons (full line) and the sum (dashed 
line) . The radial distance is scaled by the dilution parameter 
a. (b) Contour plots of the neutron perturbed density. From 
the surface to the center, the contour lines correspond to the 
density values: 0.0037, 0.0075, 0.015, and 0.023 /m" 3 . For 
larger density, colors show small density variations (3%). (c) 
Same for protons 

We observe that proton oscillations are mostly located 
at the surface of the system, which is a way to mini- 
mize the Coulomb repulsion energy. At the same time 
neutrons try to follow protons, however this motion is 
strongly affected by the isospin initial asymmetry and 
the difference between the neutron and proton orbitals 
at the Fermi energy. In fact, in neutron rich systems, 
much larger proton oscillations are observed, that is a 
way to form more symmetric fragments and hence to re- 
duce the symmetry energy. The last effect is particularly 
evident in the neutron rich 132 Sn for which the neutrons 
are more difficult to put in motion. We observe a quite 
complex structure of the unstable modes: Volume and 
surface instabilities are generally coupled and cannot be 
easily disantangled, as well as isoscalar and isovector ex- 
citations, since protons and neutrons do not move in the 
same way. Figure 3 shows the correlations between the 



amplitude of proton fluctuations p+ (r) and neutron fluc- 
tuations Pn(r), obtained at a radial distance r ranging 
from to the system radius R, for the octupole mode. 
The different panels correspond to the three Sn isotopes 
considered and to three dilutions, a, indicated in the fig- 
ure. The correlations describe closed paths, since Pp n (r) 
starts from zero at r = and goes to zero again at large 
distances. 
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FIG. 3. Correlations between the amplitude of proton and 
neutron fluctuations, for several dilutions. The full circles 
correspond to a radial distance < r < 2.4 a fm, the open 
circles correspond to 2.4 a fm < r < R. The dashed line 
indicates p+ = p+, the full line corresponds to pp = p+ Z/N. 
Results are shown for 100 Sn (a), 120 Sn (b) and 132 Sn (c). 

When protons and neutrons fluctuate whithout chang- 
ing the local chemical ratio one should have pp(r) = 
p+(r) Z/N (full line). An indentical motion of the two 
fluids should follow the diagonal Pp(r) — Pn( r ) (dashed 
line) . In N = Z nuclei (part a) the motion of protons and 
neutrons should be the same (isoscalar excitations), how- 
ever the Coulomb force introduces small differences. As 
seen in figure 2, in diluted systems (a > 1), we observe 
that, in the vicinity of the nuclear surface, density fluctu- 
ations are larger for protons than for neutrons, expecially 
in the case of neutron-rich unstable isotopes, 120 Sn and 
132 Sn (see panels (b) and (c) in figure 3), leading to a 
reduction of the asymmetry of the fragments produced 
at the surface. In the interior we observe, for these iso- 
topes, an interesting evolution from Pp{r) < p^(r) Z/N 
for a = 1 (stable modes) to pp(r) > p^i?) Z/N for the 
dilute systems. Hence we observe a change of the local 
chemical ratio, leading to a reduction of the symmetry 
energy, also in the interior of the system. This effect 
may be related to the isospin fractionation that occurs 
in unstable asymmetric nuclear matter [Ell. The pro- 
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ton migration towards the large density domains is more 
frequent than the neutron migration, which may lead to 
formation of more symmetric fragments. 

We, also, carry out calculations at finite temperature 
and determine the dilutions at which different unstable 
modes begin to appear. This allows us to specify the bor- 
der of the instability region in the density-temperature 
plane for different unstable modes. Figure 4 shows phase 
diagrams for octupole instabilities in Sn and Ca isotopes. 




00 0.05 o.io 



0.00 0.05 o.io 



0.00 0.05 0.10 0.15 



Density p(fm ) 



FIG. 4. Border of the instability region (full fine) associ- 
ated with L — 3, for Ca and Sn isotopes. The dashed line 
connects the points having the instability growth time r = 
100 fm/c. The dots are associated with r = 50 fm/c. 

Here, for simplicity, we define the density as p = po/a 3 . 
The full line indicates the border of the instability region. 
The dashed line connects points that are associated with 
the instability growth time r = 100 fm/c, and the dots 
corresponds to situations (if any) with a shorter growth 
time r = 50 fm/c. The instability region appears quite 
reduced as compared to that of nuclear matter. The lim- 
iting temperature for instability to occur is around 6 MeV 
for Sn and 4.5 MeV for Ca while it is about 16 MeV in 
symmetric nuclear matter. Heavier systems have larger 
instability region than the lighter ones. Moreover, more 
asymmetric systems are less unstable. In fact, in spite of 
the larger mass, 132 Sn is less unstable than 120 Sn. Larger 
instability growth times are obtained for 132 Sn. As seen 
from phase diagrams for 36 Ca and 40 Ca, the region of 
instability is reduced also in proton-rich systems. This 
behavior is in agreement with nuclear matter calcula- 
tions, which indicates that the instability region shrinks 
in asymmetric nuclear matter 

In conclusion, we have presented a study of the early 
development of spinodal instabilities in finite nuclear sys- 
tems by employing a quantal RPA approach. Results are 
relevant for multifragmentation studies, in fact dominant 
unstable modes determine the onset of the subsequent 
fragmentation of the system, that has to be followed us- 
ing a full dynamical calculation. We have investigated 



isospin effects on instabilities by carrying out calcula- 
tions for Ca and Sn isotopes. We find that instabilities 
are mostly of isoscalar nature, but contain also an impor- 
tant isovector component in asymmetric matter. Hence 
the liquid-gas separation of asymmetric systems is al- 
ways linked to a chemical separation inducing a fractional 
distillation of the system. The degree of instability in 
neutron-rich and also proton-rich nuclei appears reduced 
as compared to that of symmetric nuclei of comparable 
size. The instabilities are also reduced in small nuclei. Fi- 
nally we have stressed important quantum effects, such 
as the octupolar nature of the most important instability. 
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